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In high energy pp collisions J ftp's are produced di- 
rectly, from decays of higher mass charmonium states 
[X and ip(2S)], and from b quark decays. Existing exper- 
imental results in the central rapidity region from UA1 
§ at V« = 0.63 TeV, and from CDF @ and D0 § at 
•y/s = 1.8 TeV demonstrate that then measured inclu- 
sive J/ip transverse momentum distribution cannot be 
described solely by contributions from b quark decays 
and prompt production predicted by the color singlet 
model 0. In the color singlet model the charmonium 
meson retains the quantum numbers of the produced cc 
pair and thus each J/ip state can only be directly pro- 
duced via the corresponding hard scattering color singlet 
subprocess. The model predicts direct J/ip and ip(2S) 
production rates fifty times smaller than those observed 
by CDF To explain this discrepancy, a color octet 
model was introduced ^ Q . The color octet mechanism 
extends the color singlet approach by taking into account 
the production of cc pairs in a color octet configuration 
accompanied by a gluon. The color octet state evolves 
into a color singlet state via emission of a soft gluon. 
The parameters of the model were derived from a fit to 
CDF data for direct J/tp and ip(2S) production at cen- 
tral rapidity. In this article we utilize the large rapidity 
coverage of the D0 muon system to study the process 
pp —>■ J/ip + X — > p + [i~ + X in previously unexplored 
kinematical regions of small J/ip transverse momenta and 
large rapidities. We compare our results with theoretical 
predictions extended into this kinematic domain. 

The D0 detector || consists of three main systems: 
central and forward drift chambers, used to identify 
charged tracks for pseudorapidity \rj\ < 3.2; the uranium- 
liquid argon calorimeter with nearly hermetic coverage 
for \r]\ < 4; and the muon system. The detector com- 
ponent most relevant to this analysis is the Small Angle 
MUon Spectrometer (SAMUS) |9pc|] consisting of mag- 
netized iron toroids and drift tube stations on each side 
of the interaction region with pseudorapidity coverage of 
2.2 < |?7 M | <3.3 for a single muon. 

The SAMUS stations, three in each arm, consist of 
three planes of 29 mm diameter drift tubes: vertical, hor- 
izontal, and inclined at 45°. The list of tubes containing 
hits is sent to the trigger system and drift times are used 
for offline track reconstruction. Muon track reconstruc- 
tion is based on a Kalman fit jyj to the three-dimensional 
coordinates of muons passing through the SAMUS sta- 
tions, one before the toroidal magnet and two after, and 
the coordinates of the event vertex. The muon momen- 
tum resolution a p /p is about 20%, limited by SAMUS 
coordinate resolution and by Coulomb scattering in the 
calorimeter and muon toroid. 

The data were collected using the multilevel trigger 
system. The Level trigger |l^] is used to select sin- 
gle interaction events with hits in scintillator hodoscopes 
situated on both sides of the interaction region. At the 
Level 1 [O, signals from individual SAMUS tubes are 



OR'ed to provide 12 cm wide hodoscopic elements. The 
trigger requires a pattern of hits in the trigger elements 
consistent with at least one muon with transverse mo- 
mentum p 1 ^, > 3 GeV/c coming from the interaction re- 
gion. Due to high tube occupancy (~ 4%) by soft elec- 
trons and positrons, we implement a "multiplicity cut" 
at the Level 1 trigger. This cut rejects an event if the 
number of hit trigger elements in a vertically oriented 
tube plane exceeds a fixed threshold. 

The logic of the Level 1.5 trigger is similar to that of 
the Level 1 trigger, but is based on better spatial segmen- 
tation (1.5 cm vs. 12 cm). Events which pass the Level 
1.5 trigger are digitized and sent to the Level 2 software 
trigger implemented on a farm of VAX stations, where 
reconstruction of muon tracks without using drift times 
is performed. The calorimeter information is used in the 
Level 2 trigger to confirm the muon through its energy 
deposition. 

Even with the multiplicity cut, the counting rates of 
the Level 1 and Level 1.5 triggers are high in comparison 
with the allocated trigger bandwidth. To further reduce 
counting rates, we use prescales up to 10 for the dimuon 
trigger. 

In the offline analysis, we select events with one inter- 
action vertex, a single muon or dimuon trigger, and at 
least two reconstructed muon tracks. Each muon candi- 
date is required to have at least 15 hits on a track out 
of an average of 18. The energy deposition in the cells 
of the hadronic calorimeter along the muon track is re- 
quired to exceed 1.5 GeV, and to be spread contiguously 
among all five calorimeter layers. To ensure a good mo- 
mentum measurement, we require p M < 150 GeV/c and a 
minimum traverse magnetic field integral of 1.2 T-m. In 
total, 1779 events with opposite sign muon pairs and 281 
events with same sign muon pairs are selected from the 
data sample with integrated luminosity of 9.8 ± 0.5 pb _1 
|]l4| . The estimated fraction of background tracks from 
accidental hit combinations in the final data sample is 
below 1%. 

The opposite sign dimuon invariant mass distribu- 
tion M MA1 for events with transverse momentum in the 
range 1.0 < pt^P 1 < 16 GeV/c and pseudorapidity 2.5 
< < 3.7 is shown in Fig. 0. In addition to the J/ip 
signal, other contributions to the dimuon spectrum with 
Mf_if_i < 9 GeV/c 2 are expected to come from bb and cc 
production (jointly denoted as qq) with the heavy quarks 
decaying semileptonically or via sequential semileptonic 
decays, Drell-Yan production (DY), decay of light mesons 
(e.g. p, 4>, rf), and ir or K decays. 

To estimate the background and simulate the J/ip de- 
tection efficiency, we use a sample of Monte Carlo (MC) 
events from the PYTHIA 5.7 pj| and JETSET |l6| MC 
generators for each of the dimuon processes mentioned 
above, except ir and K decays. The J/ip events are gen- 
erated (assuming no J/ip polarization) with p^^ from 
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Dimuon invariant mass (GeV/c ) 

FIG. 1. The invariant mass spectrum of opposite sign 
dimuons with 1.0 < p 1 ^ < 16 GeV/c and 2.5 < |j7 mm | < 3.7. 
The hatched area indicates the J/tp signal above the sum of 
the backgrounds. 



1 GeV/c to 20 GeV/c and |?7 J /^| between 2.0 and 4.0. 
Generated dimuon events are simulated using D0 GEANT 
|l7f and mixed with minimum bias events from the data 
to simulate the combinatoric background. These simu- 
lated dimuon events are then subjected to a full trigger 
simulation and processed with the standard D0 recon- 
struction program. 

Based on MC studies we approximate the J/tp signal 
by a Gaussian function of l/M^ to account for lim- 
ited muon momentum resolution. The mass spectrum 
in Fig. [l] is fit by the sum of the J/tp signal (with the 
width and mean value as free parameters) and MC mass 
distributions for background processes (with free normal- 
ization). The number of events due to 7r and K de- 
cays is estimated from the data using like-sign dimuon 
events. The fit yields 691 ±41 J/tp events with mean 
mass (M^^) = 3.03 ± 0.03 GeV/c 2 , and standard devia- 
tion <j M = 0.56 ± 0.03 GeV/c 2 . 

The dimuon mass resolution does not allow a clear sep- 
aration of the J/tp and tp(2S) states. The fit of the in- 
variant mass distribution yields a 90% C.L. upper limit 
of the tp(2S) fraction in the signal associated with the 
J/tp of 15%. A direct measurement of inclusive tp(2S) 
production for \r]^\ < 0.6 performed by CDF shows 
that the J/tp differential cross section is approximately 
13 times larger than that of the tp(2S). 

The inclusive differential cross section of J/tp produc- 
tion is calculated from: 
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FIG. 2. The pr dependence of the J/tp differential cross 
section and its theoretical predictions (upper figure). Only 
the statistical errors are shown. The lower figure presents 
systematic uncertainties; the solid curves are the sum of all 
systematic errors, the dashed curves represent the uncertainty 
band due to J/tp polarization. The upper (lower) dashed 
curve corresponds to 100% transverse (longitudinal) polariza- 
tion. 



where L is the total integrated luminosity, ey is the J / tp 
detection efficiency, and A/y is the number of J/tp events 
in the Ap l T , Arf interval. 

To calculate the number of J/tp events, the fit to the 
mass spectrum is performed in five rf" 11 and nine p%r in- 
tervals. To reduce the errors of the fit in the high pPfP 
bins, the p^P dependence of the fraction of events at- 
tributed to J/tp is fit to a linear function and the results 
of this fit are used to obtain the number of J/tp events. 

The efficiency of J/tp detection includes acceptance, 
trigger efficiency, reconstruction efficiency, and offline 
cuts and is given by 

_ N(p l T ,Tf>) ■ e cor 
£lJ ~ N tot (p^ V i) ' 

where N [PtiV'') is the number of events in a given p z T , 
rf bin which passed all selection criteria, N tot (p^,^) 
is the total number of generated events in a bin, and 
e cor is the correction factor for effects not simulated in 
MC. The e C or includes efficiencies for the Level 2 trigger 
calorimeter confirmation of (91 ± 2)% for the dimuon 
trigger and (95 ± 1)% for the single muon trigger, and 
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TABLE I. J/ip inclusive differential cross sections Br(J/tp — > /i + fi )d 2 a/dpTdr/ (nb/GeV/c). 



Pt /iP (GeV/c) 


J fib 


2.5 - 3.7 


2.65 


2.95 


3.25 


3.55 


1.5 




187 ± 34 


— 


137 ± 55 


183 ± 57 


130 ± 45 


2.5 




77 ± 8 


— 


69 ± 18 


69.9 ± 8.1 


45.7 ± 5.3 


3.5 




21.3 ± 2.1 


21.4 ± 6.4 


23.8 ± 2.8 


17.6 ± 2.0 


15.6 ± 1.7 


4.0 




7 i A _j_ n 77 


7 _|_ 1 7 


o.UD zt U.oo 


o.oy zt u.oo 


^ /in -i- n 

D.4U zt U.oo 


5.5 




3.03 ± 0.36 


4.20 ± 0.65 


3.22 ± 0.36 


2.44 ± 0.46 




7.0 




0.667 ± 0.075 


1.07 ± 0.13 


0.758 ± 0.076 


0.77 ± 0.16 




9.0 




0.126 ± 0.023 


0.212 ± 0.033 


0.132 ± 0.021 


0.131 ± 0.039 




11.0 




0.037 ± 0.013 


0.086 ± 0.020 


0.050 ± 0.014 






14.0 




0.0057 ± 0.0023 


0.0206 ± 0.0039 


0.0064 ± 0.0020 







for offline cuts not simulated by the MC of (79 ± 4)%. 
Efficiencies for those cuts are obtained from the data and 
include (88 ± 1)% for the single vertex cut, (94 ± 3)% for 
the energy deposition cut, and (96 ± 2)% for the cut on 
the number of hits on a track. 

The measured J/ip spectrum is unfolded to correct for 
the momentum and pseudorapidity smearing using the 
technique of Ref. jl8| . The correction factors vary from 
1.7 at low p^ to 0.4 for p^ > 8 GeV/c. 

The calculated differential cross section is fit to an ex- 
ponential function. The results are used for interpola- 
tion of the cross sections from average values of (p 1 t) 
and (l^y-' l) to the centers of the selected intervals. The 
inclusive differential J/ip cross section averaged over a 
rapidity range of 2.5 < \r) J l^\ < 3.7 is shown in Fig. |. 
Results for finer rapidity bins are collected in Table 1. 
The uncertainties quoted there are statistical only. 
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10 r * D0 (Phys. Lett. B370, 239(1996)) 
; o CDF (Phys. Rev. Lett. 79, 572(1997)) 
— Prompt (color octet+color singlet) + b 



TABLE 11. Systematic errors of the J/ip cross sections. 



10 



0.5 



1.5 



2.5 



Source 



Systematic Error 



Unfolding procedure 


15% 


J/ip background determination 


7.2% - 30% 


J/ip detection efficiency 


7% 


Level 1 multiplicity cut 


6% 


Total integrated luminosity 


5.4% 


ip(2S) contamination 
d 2 a/dpTdr/ total 


+0% 
-5% 


20% - 36% 



Averaging over 2.5 < 77 J/l/ ' < 3.7 


596 - 


30% 


da/dprAr/ total 


21% 


-47% 



The largest (> 2%) systematic uncertainties are sum- 
marized in Table |J. The contribution from the unfolding 
is derived from comparison with the bin-by-bin unfold- 
ing technique |is|| . The uncertainties in the determina- 
tion of the background and averaging cross section over 
the SAMUS pseudorapidity acceptance vary for different 
pij,^ and are caused by uncertainties in the parameter- 
ization of the data. The difference in the parameters of 
measured and generated J/ip mass distributions as well 
as the accuracy of the spectrometer description in the 
detector simulation are used to estimate the J/ip detec- 
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FIG. 3. The pseudorapidity dependence of the J/ip pro- 
duction cross section with pr > 5 GeV/c (upper points and 
curve) and pr > 8 GeV/c (lower points and curve). The 
error bars are statistical and systematic errors (polarization 
uncertainties not included) summed in quadrature. 



tion efficiency uncertainty. The uncertainty due to the 
Level 1 multiplicity cut was determined by varying the 
threshold of this cut by one trigger element. The results 
in Table | are obtained for the case of zero J/ip polar- 
ization. The additional uncertainty up to 1^5% due to 
possible J/ip polarization is shown in Fig. ^| along with 
the p 3 J^ dependence of the total systematic error. 

In Fig. [2] we compare the J / ip cross section with cur- 
rent models of charmonium production. For J/ip from 
b quarks we use the NLO QCD predictions [^o) with the 

renormalization/factorization scale /1 = + p^ 2 , 

where mi and p^ are the parent b quark mass and trans- 
verse momentum, respectively. The scale is chosen to 
match theory predictions to the published D0 b quark 
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cross sections in the central rapidity region ||. We use 
ISAJET 1 21 1 to fragment b quarks into J/tp. The color 
octet and color singlet contributions to the direct J/tp 
production and radiative x decays are taken from Ref . . 
The term representing the direct J/tp production is in- 
creased by 12% to account for the contribution from 
ip(2S) decays @. 

Fig. H shows the pseudorapidity dependence of the 
measured J/tp cross section for p 3 J^ > 5 and 8 GeV/c 
along with the corresponding central rapidity measure- 
ments of D0 ||] and CDF Within uncertainties, 
the color octet model plus b quark decays describe the 
r\ dependence of the inclusive J/tp production in the full 
rapidity region. 

In conclusion, we have made the first measurement of 
inclusive J/tp production in the forward rapidity region 
2.5 < \j] J ^\ < 3.7 in pp collisions at ^ = 1.8 TeV. The 
data show good agreement with the theoretical predic- 
tions based on b quark decays and the color octet model 
of direct charmonium production. 
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